Mott transition in one dimension: Benchmarking dynamical cluster approaches 
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The variational cluster approach (VCA) is applied to the one-dimensional Hubbard model at zero 
temperature using clusters (chains) of up to ten sites with full diagonalization and the Lanczos 
method as cluster solver. Within the framework of the self-energy- functional theory (SFT), dif- 
ferent cluster reference systems with and without bath degrees of freedom, in different topologies 
and with different sets of variational parameters are considered. Static and one-particle dynamical 
quantities are calculated for half-filling as a function of U as well as for fixed U as a function of 
the chemical potential to study the interaction- and filling-dependent metal-insulator (Mott) transi- 
tion. The recently developed Q-matrix technique is used to compute the SFT grand potential. For 
benchmarking purposes we compare the VCA results with exact results available from the Bethe 
ansatz, with essentially exact dynamical DMRG data, with (cellular) dynamical mean-field theory 
and full diagonalization of isolated Hubbard chains. Several issues are discussed including conver- 
gence of the results with cluster size, the ability of cluster approaches to access the critical regime of 
the Mott transition, efficiency in the optimization of correlated-site vs. bath-site parameters and of 
multi-dimensional parameter optimization. We also study the role of bath sites for the description 
of excitation properties and as charge reservoirs for the description of filling dependencies. The 
VCA turns out to be a computationally cheap method which is competitive with established cluster 
approaches. 

PACS numbers: 71.10.-w, 71.30.+h, 71.10.Fd 



I. INTRODUCTION 

One of the most fascinating effects of strong inter- 
actions among itinerant electrons is insulating behav- 
ior that is induced by electron correlations An un- 
derstanding of the Mott insulating state and also of the 
transition from a correlated metal to a Mott insulator is 
important for many transition-metal oxides including the 
parent compounds of cuprate-based high-temperature 
superconductors.— The Mott transition is generically de- 
scribed using lattice models with purely local interac- 
tions, such as the single-band Hubbard model 4^ Cor- 
relations, however, are generally non-local or even long- 
ranged. It is a highly non-trivial question to what ex- 
tent the Mott transition is dominated by local correla- 
tions and whether or not local approximations, i.e. ap- 
proximations neglecting non-local correlations, are able 
to capture the essence of the Mott physics. 

The implications of the dynamical mean-field 
theory^&ii as a distinguished local approximation, have 
been worked out in detail in the past and have been com- 
pared with numerous experiments on transition-metal 
oxides. The DMFT is a mean-field approach in the sense 
that the effects of non-local spin or charge two-particle 
correlations on the one-particle excitation spectrum are 
neglected. In case of layered, essentially two-dimensional 
compounds, for example, this is probably a too strong 
approximation. 

Different generalizations of DMFT have been sug- 
gested in the past to overcome this limitation. 10 ' 11 ! 12 ! 13 
Cluster extensions^ which restore the effects of non-local 



correlations step by step with increasing cluster size, are 
conceptually simple and interesting approaches in this re- 
spect. The idea is that, besides the local (temporal) cor- 
relations, it is the effect of the short-range correlations, 
treated exactly in a cluster approach, which dominates 
the physics of the Mott transition or at least the physics 
of the Mott-insulating state. As the importance of non- 
local correlations is expected to increase with decreasing 
lattice dimension, the most stringent test for a cluster ap- 
proach consists in its application to the one-dimensional 
lattice. Ideally, starting from a (dynamical) mean-field 
picture of the transition, the decisive step forward should 
be done with the smallest clusters already, while a fur- 
ther increase of the cluster size should add qualitatively 
unimportant corrections only. Since exact results from 
the Bethe ansatz are available for the Hubbard model in 
one dimension^ this model can very well be used to test 
this idea and to benchmark different cluster approxima- 
tions. 

Apart from true extensions of DMFT^ i 16 ' 17 ^ where 
a small cluster with L c correlated sites and a continuum 
(n s = oo) of uncorrelated ("bath") sites is considered, 
there are also dynamical cluster approaches without any 
bath degrees of freedom (i.e. n s = 1), such as the cluster- 
perturbation theory (CPT) 19 i 20 and the variational clus- 
ter approach (VCA) i2i The self-energy- functional theory 
(SFT) ) 21 ' 22 provides a unified framework for all these dif- 
ferent cluster approximations which are characterized by 
different L c and n s . Therefore, the purpose of this pa- 
per is to apply (within the SFT) different cluster ap- 
proaches to the one-dimensional Hubbard model and to 
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study the interaction- and the filling-controlled transition 
for benchmarking purposes. 

An obvious question concerns the rate of convergence 
to the exact solution which is approached in the L c = oo 
limit. Depending on the cluster scheme used and on the 
quantity of interest, an exponential or power-law depen- 
dence on L c can be inferredi 23 i 24 ' 25 i 26 These considera- 
tions, however, apply to the largc-L c limit only while for 
practical purposes the quality of a given approximation 
for small L c is much more important and can only be 
estimated a posteriori. 

It is also unclear whether or not the inclusion of bath 
sites makes a cluster approach more efficient, i.e. speeds 
up the convergence to the exact solution (see the discus- 
sion in Ref. 1211 ). Because of the exponential growth of 
the Hilbert space with increasing n s , this is of particu- 
lar importance for methods that are based on a full (or 
Lanczos) diagonalization of the effective cluster model. 

A simple and frequently used2£ cluster approximation 
is the "direct" cluster approximation where quantities 
for the infinite system are approximated by those of a 
finite cluster without any embedding of the cluster into 
a medium that mimics the disregarded cluster environ- 
ment. Due to particle-number conservation and due to 
the finite (and usually small) cluster size L c , it is in- 
convenient to study filling dependencies using the direct 
cluster approach. The filling-controlled Mott transition, 
in particular, is hardly accessible in this way. Another 
interesting question is therefore, if and how a continuous 
filling dependence can be achieved by embedded, self- 
consistent or variational cluster approaches. 

In the following Sec. [Til we briefly list some well-known 
properties of the Hubbard model in one dimension which 
are relevant for our study. Sec. IIIII provides a brief dis- 
cussion of cluster approaches employing the exact diago- 
nalization method at zero temperature. Some details of 
the SFT and of our numerical approach are given in Sec. 
IIVI The numerical results for the Mott insulator at half- 
filling and for the filling-controlled Mott transition are 
presented and discussed in Sees. IVl and I VI) respectively. 
Sec. IVIII summarizes our main results. 



of the local Coulomb repulsion, and fi is the chemical 
potential. 

The ground-state energy Eq (of H + fj,N) can be cal- 
culated exactly^ 4 - by means of the Bethe ansatz. For 
fj, = U/2 the model is particle-hole symmetric and half- 
filled (n = 1). If L denotes the number of sites, we have 
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where Jo, Ji are Bessel functions. While the system is 
metallic for U — 0, a Mott-insulating state is found for 
any U > as can be seen from the exact expression 2 ^ for 
the single-particle gap: 
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The gap is finite for U > but exponentially small in 
the limit U — > 0, i.e. A ~ exp(— 1/U). Opposed to the 
dynamical mean-field scenario, there is no finite critical 
U for the Mott transition. 

The filling-controlled Mott transition can be charac- 
terized by the charge susceptibility ( "compressibility" ) 
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At U — the compressibility is proportional to the tight- 
binding density of states at the Fermi energy and is there- 
fore finite for all fillings, except for n = and n = 2 be- 
cause of the van Hove singularities at the lower and upper 
band edge. For any finite U, the compressibility must 
vanish in the Mott-insulating phase as for n = 1 the one- 
particle excitation spectrum is gapped. Approaching the 
Mott insulator from the metallic side (n — > 1), however, 
the compressibility behaves discontinuously and even di- 
verges. For n < 1 but close to half-filling, it is given 
by^ 



K = aS 



(5) 



II. THE MOTT TRANSITION IN THE 
ONE-DIMENSIONAL HUBBARD MODEL 

In one dimension, the single-band (grand-canonical) 
Hubbard model is given by 

H = -t^ (c^a-ia + H.c.J - fi^ma + U^ n^riii . 

i<7 icr icr 

(1) 

Here, C{ a annihilates an electron at the site i with spin 
projection a =T,|- Furthermore, = c ia Ci a is the 
occupation-number operator, the ground-state average of 
which is the site- and spin-independent filling n = (n^). 
We consider nearest-neighbor hopping only and set t = 1 
to fix the energy scale. Finally, U denotes the strength 



where 8 = 1 — n is the hole concentration and a > a In- 
dependent constant. This implies that close to half-filling 
<5, as a function of the chemical potential, approaches 
(5 = with an infinite slope: 

5 oc V^(O) - W) ■ (6) 

Here, we note that the Mott physics in dimension D = 
1 stands in marked contrast to the DMFT (or D = oo) 
scenario^ DMFT predicts a Mott-insulating state with a 
vanishing k only for interaction strengths U larger than 
a finite critical value U c . Furthermore, for U > U c the 
compressibility is found to stay finite, < k < oo, when 
approaching the Mott insulator from the metallic state 
off half-filling as n — > 1. 
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III. CLUSTER APPROACHES USING THE 
LANCZOS TECHNIQUE 

It is obviously interesting whether or not this qualita- 
tive difference in the physics of the Mott transition can be 
captured by means of a cluster extension of the DMFT. 
This question has been tackled recently with the help of 
cellular DMFT (and with the periodized C-DMFT) for 
half-filling by Bolech et al^ and for the filling-dependent 
transition by Capone et al.— For studies at zero temper- 
ature the Lanczos technique^ is a powerful method to 
treat the effective cluster problem within C-DMFT and 
was also employed in Refs. I30ll3ll The effective cluster 
Hamiltonian is given by 

H ' = -tEE(cU-iaH,.) 

i—2 cr 

i=l a i=l u 

n,(i) 

+ [ £ i^ a \ Vi a a ^ + ( V iVi C L a Wi° + H - C -) 



v i= 2 



(7) 



Here n s (i) — 1 is the number of uncorrelated sites per cor- 
related site i. n s (i) may vary along the chain. Due to the 
exponential growth of the Hilbert space dimension with 
the system size, calculations are basically limited to a to- 
tal number of J2t=i n s(i) ~ 10 s ^es only, if L c denotes 
the number of correlated sites (the cluster "size"). Note 
that the reference system has to be solved repeatedly to 
find a stationary point or to achieve self-consistency, re- 
spectively, and that due to open boundary conditions a 
few general symmetries can be exploited only. For a small 
cluster with L c = 4, for example, this implies a limita- 
tion to less than n s = 3 local degrees of freedom, i.e. less 
than two bath sites per correlated site, if n s (i) = n s is 
taken to be constant as usual. While a true solution of 
the C-DMFT self-consistency equation actually requires 
a continuum of bath sites (at each correlated site), i.e. 
n s =oo, the convergence with respect to n s is expected^ 3 - 
to be exponentially fast. This makes calculations with 
small n s feasible. 

However, there are two conceptual drawbacks of the 
exact-diagonalization (Lanczos) approach to C-DMFT: 
(i) Clearly, the determination of the bath parameters 
Vi V is of crucial importance for small n s (and small L c ). 
One possible prescription is to fix the parameters by min- 
imization of a suitably defined distance between the hy- 
bridization function of Eq. ([7|) and the one given by the 
self-consistency equation. The choice of the quantity that 
is "projected" as well as the distance measure, however, 
are more or less ad hoc and by no means unique, (ii) 
Within the C-DMFT the one-particle energies of and the 
hopping between the correlated sites are fixed by their 
values in the original Hubbard model. This may be seen 



as a limited flexibility for the determination of the (in a 
certain sense) optimal effective cluster model. 

The variational cluster approximation (VCA)^! or, 
more generally, the self-energy-functional theory22 does 
not suffer from these shortcomings: (i) The bath param- 
eters of the effective cluster model (the "reference sys- 
tem" ) are fixed in a unique way by demanding the grand 
potential of the system to be stationary with respect to 
those variations of the self-energy that are induced by 
varying the bath parameters. This prescription is dis- 
tinguished by the fact that it ensures thermodynamical 
consistency of the results; 34 i 35 All approximate quanti- 
ties of the theory derive from an approximate but ex- 
plicitly given thermodynamical potential. Opposed to 
C-DMFT/ED this consistency is achieved for any L c , n s 
and not only in the continuum limit n s — > oo. (ii) There 
is more flexibility in the choice of the reference system: 
Within the SFT it is possible to vary all one-particle pa- 
rameters of the reference system including those referring 
to the original correlated sites. Furthermore, one is by no 
means forced to attach a bath to each of the correlated 
sites. A physically motivated choice is to consider bath 
sites at the cluster boundaries only, for example. 



VARIATIONAL CLUSTER APPROACH 
USING Q MATRICES 



IV. 



The SFT is described in Refs. |2Tj|22|3J|35||36||37|. The 
main idea is to express the grand potential of the original 
model as a functional of the self-energy, fl = f2[X], such 
that the exact self-energy is given as a stationary point, 
<5f2/<ffi = 0. Trial self-energies arc taken from a refer- 
ence system with the same (Hubbard) interaction but 
with a modified one-particle part. If the Hamiltonian of 
the original system, H = H (t) + Hi(U), consists of a 
free part with parameters t and an interaction term with 
parameters U, the most general Hamiltonian of the refer- 
ence system has the form H' = Ho(t') + Hi(U) with arbi- 
trary t' . Fig. [T] shows the original one-dimensional Hub- 
bard model with nearest-neighbor hopping t as well as 
various reference systems considered for our calculations. 
The trial self-energy is parametrized by the set of one- 
particle parameters of the reference system: £ = E(i'), 
and variations of the trial self-energy are considered that 
are due to variations of t' . i.e. one has to solve: 



d»[S(t')] 
di' 







(8) 



The decisive point is that O [£(£')] can be evaluated ex- 
actly for reference systems that allow for a (numerically) 
exact computation of the single-particle Green's function. 
In case of a finite (small) cluster or chain and a finite 
(small) number of bath sites, this can be achieved by full 
diagonalization or with the help of the Lanczos method.— 
All what is needed in a practical calculation is the one- 
particle Green's function of the reference system. If L c > 
1, i.e. in case of the variational cluster approximation 
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(VCA), this is the Green's function of a set of decoupled 
clusters. The Green's function for a single cluster, 

u> — 10 

is given in terms of poles aj' m and corresponding weights 
QamQmp- The poles and the Q-matrices^i can be read 
off from the standard Lehmann representation^ Note 
that Q is a non-quadratic matrix: a = (i, a) refers to a 
one-particle orbital of the cluster while m = (r, s) refers 
to a single-particle excitation between two eigenstates |s) 
and |r) of the cluster Hamiltonian H' with excitation 
energy w' m =E' r - E' s . We have QQ f = 1 ^ Q^Q. 
For T = the SFT grand potential is then given by: 

Q[s(t')\ = n'+J2u m e(-co m )-J2"' m ®(-0- (io) 

m m 

Here O' is the grand potential of the reference system, 
0(w) is the Heaviside step function, and u) m are the 
poles of the (VCA approximation for the) one-particle 
Green's function of the original system. They can easily 
be obtained^! as the eigenvalues of the matrix 

M = A + QWQ , (11) 

with A mn = uj' m 8 mn and V = t — t'. Typically, V in- 
cludes the inter-cluster hopping, shifts of one-particle en- 
ergies and, in the case of bath sites, further hybridization 
terms. 

If a reference system with bath sites is considered, it 
is convenient to formally include these bath sites also in 
the original system where they are, however, completely 
decoupled from the correlated sites such that physical 
quantities remain unchanged. This has the advantage 
that t and t! have the same matrix dimension, and also 
the Hamiltonians of the original and of the reference sys- 
tem, H and H', operate on the same Hilbert space. 

The dimension of M is given by the number of poles 
of G with non- vanishing spectral weight. If this number 
is not too large, the Q-matrix technique is a very simple 
means for the evaluation of the self-energy functional. 
As there is no frequency integration involved, neither on 
the real axis where an additional broadening parameter 
must be used, nor on the imaginary axis where a high- 
frequency cutoff must be introduced and the remaining 
tail must be controlled, the method is also very accurate. 

For larger clusters not accessible to full diagonaliza- 
tion, we employ the band Lanczos method to compute A 
and Q?^ This variant of the Lanczos algorithm ensures 
that different elements G' a/3 (ui) have the same set of poles, 
i.e. the same u' m independent of a, (3. The dimension of 
the matrix M is then given by the number of iteration 
steps in the Lanczos procedure. As the results usually 
converge very fast, it is sufficient to consider about 100 
steps only. This is regularly checked in our calculations. 
For small clusters the results have been compared with 
those obtained by full diagonalization and found to agree 
within numerical accuracy. 



For a given reference system one should in principle 
vary all one-particle parameters t' to get the optimal re- 
sult. It is much more convenient, however, to restrict 
oneself to a small number of physically motivated param- 
eters to be optimized. This avoids complications arising 
from a search for a stationary point in a high-dimensional 
parameter space. In most cases, as will also be demon- 
strated below, it is fully sufficient to consider a few vari- 
ational parameters only which are suggested by the ge- 
ometry of the reference system in an obvious way. The 
reference systems considered here as well as the corre- 
sponding variational parameters taken into account are 
shown in Fig. Q] 

The system of interest is the D = 1 Hubbard model. 
However, for practical purposes it is more convenient to 
consider a Hubbard chain consisting of a finite number of 
sites L with periodic boundary conditions as our original 
system. For the actual calculations we used L = 1000 — 
2000 sites. This is fully sufficient to ensure that all results 
shown below are independent of L. 

Stationary points are determined using different nu- 
merical strategies^ One-dimensional parameter opti- 
mization is performed by iterative bracketing of maxima 
and minima. For more than one variational parameter, 
the SFT grand potential is usually not extremal but has a 
saddle point. Given a certain characteristic of the saddle 
point, this can be found by iterated one-dimensional opti- 
mizations - a strategy that has been found to be useful for 
two or three parameters. In case of higher-dimensional 
parameter spaces, the downhill simplex method is applied 
to find local minima of \d£l[£(t')]/dt'\ 2 from which (if 
there is more than one) only those are retained for which 
£l[S(t')] has a vanishing gradient. For most situations^ 
the minimal grand potential distinguishes the thcrmody- 
namically stable phase if there is more than a single sta- 
tionary point. In all examples discussed below, however, 
this has not been an issue or turned out to be straight- 
forward. 



V. MOTT-INSULATING PHASE FOR 
HALF-FILLING 

One of the advantages of the variational cluster ap- 
proximation (and of the SFT in general) consists in its 
flexibility to construct approximations of different qual- 
ity and complexity. The most simple case is given by a 
reference system consisting of decoupled clusters (chains) 
with L c correlated sites each where only the intra-clustcr 
nearest-neighbor hopping t' is considered as a variational 
parameter, see Fig. [1] A. The hopping is assumed to be 
the same for all pairs of nearest neighbors. This im- 
plies that the parameter space is one-dimensional, and 
an extremumoi Sl(i') = £![X(£')] defines the physical self- 
energy within this approximation. Clusters with an even 
number of sites L c are preferred to avoid a Kramers- 
degenerate ground state (for half-filling) and odd-even 
effects when comparing results for different L c with each 
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FIG. 1: (Color online) The original system (D = 1 Hubbard 
model, n.n. hopping) and different reference systems consid- 
ered in this study together with the corresponding variational 
parameters being optimized. See text for discussion. 
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FIG. 2: (Color online) [/-dependence of the optimal intra- 
cluster hopping t' for different cluster sizes L c as indicated. 
VCA calculations for (j, = U/2 (half-filling) using the reference 
system displayed in Fig.[T] A. The physical n.n. hopping t = 1 
sets the energy scale. Inset: Optimal t' as function of U 2 for 
U —>0. 



A. Variation of hopping parameters 



other. 

Another variational parameter suggesting itself is the 
on-site energy. A homogeneous shift of all infra-cluster 
on-site energies acts like a separate cluster chemical po- 
tential. As has been discussed in Ref. HH, this is actually 
one of the most important parameters as its optimization 
guarantees thermodynamical consistency with respect to 
the particle number. Here, we start our discussion with 
the Mott- insulating phase at half- filling (/x = U/2). In 
this case, the optimal value of the cluster on-site ener- 
gies is already predicted by particle-hole symmetry and 
is given by the on-site energy of the sites in the original 
model (which is set to zero). This can also be reproduced 
explicitly within the VCA: It turns out that f2[E(£')] is 
always at a maximum for vanishing on-site energies. 

We also ignore fictitious symmetry-breaking fields cou- 
pling to the local spin or the local charge density. Since 
those coupling terms belong to the one-particle part 
of the Hamiltonian, the corresponding field strengths 
can in principle be considered as additional variational 
parameters^ A finite value of the optimized field would 
indicate long-range spin or charge order which, however, 
is absent in one dimension or, as concerns e.g. ferromag- 
netism, is disregarded here. As discussed in Ref. Hi], the 
absence of antifcrromagnctic order in the D = 1 Hub- 
bard model at half-filling is respected by the VCA for 
sufficiently large clusters. 



A non-trivial result, namely t' 7^ t, is found when op- 
timizing the nearest-neighbor hopping, see Fig. [2l The 
physical idea behind this approximation is that switching 
off the inter-cluster hopping, which generates the approx- 
imate self-energy, can partially be compensated for by en- 
hancing the mtra-cluster hopping. This is in fact seen in 
the figure: The optimal t' is larger than the physical hop- 
ping. The trends found for different cluster sizes L c and 
for different U corroborate this interpretation: The larger 
the cluster the smaller is the necessary compensation (see 
Fig. [2]). Furthermore, it is reasonable that in case of a 
stronger interaction and thus more localized electrons, 
switching off the inter-cluster hopping is less significant. 
The strongest approximation of the self-energy is there- 
fore generated by the smallest cluster (L c = 2) and in the 
limit U — > 0. This is indicated by a strong (more than 
100%) enhancement of t' compared to t. 

On the other hand, even a "strong" approximation for 
the self-energy (measured as a strong deviation of t' from 
t) becomes irrelevant in the weak-coupling limit because 
the self-energy must vanish for [7 = 0. It is therefore 
not surprising that the VCA exactly recovers the U = 
limit. The approximate VCA Green's function, which 
can be expressed as G(u>) = (Go(w) -1 — S(cj)) -1 in terms 
of the optimized self-energy and the free lattice Green's 
function, becomes exact for U = 0. The same holds 
for the SFT grand potential at the stationary point fi 
or for the ground-state energy Eq = ft + fi(N). The 
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FIG. 3: (Color online) 17-dependence of the VCA ground- 
state energy per site for cluster size L c = 2 and L c = 10 (see 
Fig-IU A). The exact (Bethe ansatz) result of Eq. Q is shown 
for comparison. 



latter is shown in Fig. [3] as a function of U in comparison 
with the exact (Bethe ansatz) result of Eq. @. Both, 
VCA calculations for the smallest (L c = 2) as well as for 
much larger (L c = 10) cluster size, correctly reproduce 
the U = limit while for strong interactions there are 
deviations. As expected the L c = 10 calculation provides 
a much better approximation. 

Fig.[4]dcmonstrates how the U = limit is approached. 
For strong interaction U = 10 the SFT grand potential 
f2(t') = f2[E(t')] is at a minimum for t 1 as t. Upon de- 
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FIG. 4: (Color online) SFT grand potential fi(t') = n[S(i')] 
(constantly shifted by fi{N)) per site and as a function of 
the intra-cluster nearest-neighbor hopping for L c = 10 and 
different U (fi — (7/2). Arrows indicate the respective optimal 
t'. 
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FIG. 5: (Color online) VCA ground-state energy per site for 
(7 = 4 (top) and (7 = 8 (bottom) for different cluster sizes L c 
as a function of 1/L C compared to the exact (BA) result and 
the direct cluster approach. 



creasing U , the optimal t' more and more deviates from 
the physical t. At the same time, however, the SFT grand 
potential f2(f') becomes flatter and flatter, and for U — > 
the optimal t' is completely irrelevant as S(w) = for 
any t'. 

For finite U the quality of the cluster approximation is 
determined by the cluster size L c . The dependence of the 
VCA ground-state energy Eq on L c turns out to be quite 
regular. Plotting the results for fixed U as a function 
of 1/L C allows to recover the exact ground-state energy 
by extrapolation to 1/L C = 0. This is demonstrated in 
Fig. It is worth mentioning that the VCA represents 
a considerable improvement as compared to the "direct" 
cluster approach where Eq is simply approximated by 
the ground-state energy of an isolated Hubbard chain 
(with open boundary conditions). Convergence to the 
exact result is clearly faster within the VCA. As can be 
seen by comparing the trends for U — 4 and U = 8 
in Fig. [51 this advantage is more pronounced for weaker 
interactions which is explained by the fact that, opposed 
to the VCA, the direct cluster approach is approximate 
even for (7 = 0. 

Cluster-perturbation theory (CPT)^^ can be under- 
stood as being identical with the VCA provided that the 
SFT expression for the grand potential is used and pro- 
vided that isolated clusters are used as reference system 
and no parameter optimization at all is performed. (It 
should be mentioned that this implies the inability of 
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VCA and from the direct cluster approach for L c = 2 and 
L c — 10 in comparison with the exact result of Eq. © ■ VCA 
calculations using the reference system shown in Fig. [T] A. 



CPT to describe antiferromagnetic order for D = 2 and 
T = 0, for example). As can be seen from Fig. 2J there is 
a gain in binding energy due to the optimization of t' , i.e. 
Q(t') < Cl(t) for the optimal t' . This means that the VCA 
improves on the CPT result. One should note, however, 
that on the energy scale used in Fig. for example, this 
binding-energy gain would hardly be visible. 

The VCA value for Eq, though in Fig. \5\ always higher 
than the exact result, does not represent an upper bound 
to the exact ground-state energy a prioriM- To our knowl- 
edge there is no general proof that the self-energy func- 
tional is convex or "variational" despite several recent 
efforts . 42 i 43 i 44 This must be seen as a disadvantage as 
compared to the direct cluster method which, via the 
Ritz variational principle and in the case of open bound- 
ary conditions is easily shown to provide strict upper 
bounds. However, this disadvantage appears to be in- 
herent to all hitherto known variational principles that 
are not derived from the Ritz principle. 

The VCA derives from a dynamical variational princi- 
ple based on the one-particle self-energy as the basic vari- 
able. One should therefore expect that the approach is 
able to predict one-particle excitations significantly bet- 
ter than the direct cluster method. Here, for the discus- 
sion of the Mott insulator, the focus is on the insulating 
single-particle excitation gap A, the exact [/-dependence 
of which is given by Eq. ([3]). Using the Q-matvix ap- 
proach, we get the poles of the one-particle Green's func- 
tion with finite spectral weight by diagonalization of the 
matrix M in Eq. (fTTj) . The difference between the low- 
est pole in the electron-addition part and the highest pole 
in the electron-removal part of the spectrum defines A. 
As can be seen in Fig. [6l the VCA results for L c = 2 
and L c = 10 considerably improve upon the results of 
the direct cluster method. For intermediate and strong 
couplings, the VCA gap calculated for L c = 10 is close 



to the exact result, and taking into account the L c = 2 
calculation in addition, a finite-size scaling appears to be 
possible. In the weak-coupling limit (U < 2), however, 
an increase of the cluster size apparently does no longer 
lead to a significant improvement. Although the VCA 
gap approaches A = for U — > 0, there is a clear overes- 
timation as compared to the exact result with a relative 
error that even diverges for U — > 0. 

For a more detailed discussion of the critical point 
U = consider the inset in Fig. [21 One can see that for 
U — > 0, the optimal infra-cluster hopping quadratically 
approaches a finite value: t' = t' + const, x U 2 + 0(U 3 ). 
This implies that cluster eigenenergies and thus exci- 
tation energies as well as cluster eigenstates and thus 
spectral weights depend (for L c < oo) analytically on 
U for U — ► 0. Consequently, the same holds for the 
VCA Green's function since this can be expressed in 
terms of the cluster Green's function G as G(u>) = 
(G'(w) - — (t — 1 an d since the matrix inversion 
involves finite blocks only due to the remaining super- 
lattice translational symmetry of the reference system. 
One-particle correlation functions, the ground-state en- 
ergy etc. are therefore analytical in U for U — > within 
the VCA. The same holds for the one-particle excita- 
tion gap while the exact gap is non-analytic at U = 
(cf. Eq. ©). That this non-analyticity cannot be repro- 
duced within the VCA, should be interpreted as a rather 
general failure that is inherent to any cluster concept. 
Qualitative changes at a critical point resulting from the 
limiting process L c — > oo are beyond a scheme based on 
finite clusters. 

So far the discussion has been restricted to calculations 
using a single variational parameter. More parameters 
can be useful for different reasons. First, we note that 
the optimal self-energy provided in a real-space cluster 
technique does not reflect the full translational symmetry 
of the original lattice problem and that finite-size effects 
are expected to be the most pronounced at the cluster 
boundary. This suggests to use reference systems with 
site- or bond-dependent variational parameters. For the 
case of particle-hole symmetry, obvious choices are dis- 
played in Fig.Q] B where the intra-cluster hopping at the 
edges of the chain are allowed to take a different value, 
and in Fig. [TJ C where more or all hopping parameters 
are varied independently. 

Fig. [7| shows the numerical results for [7 = 4. We find 
that the optimal hopping varies between different nearest 
neighbors within a range of less than 10%. At the chain 
edges the optimal hopping is enhanced to compensate 
the loss of itinerancy due to the switched-off inter-cluster 
hopping within the VCA. With increasing distance to the 
edges, the hopping quickly decreases. Quite generally, 
the third hopping parameter is already close to the phys- 
ical hopping t. Looking at the L c = 10 results where 
all (five) different hopping parameters have been varied 
independently (orange circles), one can see the hopping 
to slightly oscillate around the bulk value reminiscent of 
surface Fricdcl oscillations. 
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FIG. 7: (Color online) Optimized hopping parameters for ref- 
erence systems shown in Fig. [TJ A, B and C for U = 4 and 
different cluster sizes ranging from L c = 4 to L c = 10 as 
indicated (physical hopping set to t = 1). Black: hopping as- 
sumed to be uniform (A). Red: two hopping parameters var- 
ied independently, the hopping at the two cluster edges and 
the "bulk" hopping (B). Blue: hopping at the edges, next to 
the edges and bulk hopping varied. Green: four hopping pa- 
rameters varied. Orange: five hopping parameters varied (for 
L c = 10 this corresponds to C). 



The optimal SFT grand potential is found to be lower 
for the inhomogeneous cases as compared to the homo- 
geneous (black) one. Generally, the more variational pa- 
rameters are taken into account the higher is the decrease 
of the SFT grand potential at optimal parameters. How- 
ever, the binding-energy gain due to inhomogeneous hop- 
ping parameters is much smaller compared to the gain 
obtained with a larger cluster. Likewise, there is merely 
a marginal improvement as concerns the single-particle 
gap. 

Considering an additional hopping parameter i p b c link- 
ing the two chain edges as a variational parameter (Fig. 
[TJ D), always gives a minimal SFT grand potential at 
^pbc = 0. This implies that open boundary conditions 
are preferred as compared to periodic boundary condi- 
tions (which would be given by a stationary point at 
^Dbc = !)• The issue has already been discussed in Ref. 

M 

With the reference system Fig. [TJ E we can check 



whether or not a magnetic frustration develops in the ref- 
erence system. A hopping t' 2 between next-nearest neigh- 
bors leads in the Heisenberg limit U — > oo to an antifer- 
romagnetic next-nearest-neighbor exchange J2 and thus 
to a frustration of antiferromagnetic (short-range) order. 
This would partially compensate the residual mean-field 
character of the VCA with respect to magnetic proper- 
tics. At the same time, however, particle-hole symmetry 
would be violated. It turns out, however, that the SFT 
grand potential has a saddle point with t' 2 = 0. (It is at 
a minimum w.r.t. t\ and at a maximum w.r.t. t 2 ). 

A third-nearest-neighbor hopping would not lead to 
frustration and would also respect particle-hole symme- 
try. Optimization of an L c = 6-site cluster at U = 4 as 
indicated in Fig. [TJ F yields an optimal nearest-neighbor 
hopping t'i ~ 1.04 and third-nearest-neighbor hopping 
t' 2 ~ —0.02. This shows that hopping parameters that 
are not present in the original system can get a finite 
value when treated as variational parameters in the ref- 
erence system. The corresponding decrease of the SFT 
grand potential is marginal, however. Consequently, we 
disregard such variational parameters in the following. 



B. Bath degrees of freedom 

A different possibility to increase the number of varia- 
tional parameters is to introduce additional uncorrelated 
("bath") sites. As there is no Hubbard interaction on 
the bath sites, the interaction part of the Hamiltonian 
is left unchanged, as it is necessary for an allowed refer- 
ence system within the SFT*22, Note that the trial self- 
energy T,ij(uj) is still labelled by the correlated sites only. 
We consider reference systems where all or some of the 
original correlated sites are coupled to bath sites via a 
hopping ("hybridization") of strength V. For each corre- 
lated site i the different hybridization parameters Vi Vi as 
well as the one-particle energies of the bath sites Ei Vi for 
Vi = 2, n s (i) can be treated as variational parameters. 
Here n s (i) — 1 is the number of bath sites for a given cor- 
related site i. The inclusion of bath sites improves the 
description of temporal instead of spatial degrees of free- 
dom. For L c = 1 one recovers the DMFT, for L c > 1 the 
cellular DMFT in the limit n s (i) = n s — > 00^ Calcula- 
tions using the Lanczos method are feasible, however, for 
small n s only. 

Particle-hole symmetry considerably reduces the num- 
ber of variational parameters that have to be varied in- 
dependently. For a single bath site (n s (i) = 2), the 
bath on-site energy is pinned to the chemical potential, 
£ = fj, = U/2, and only the hybridization V is free. For 
n s {i) = 3 we have e = fi ± Ae with a variational param- 
eter Ae. Both bath sites couple with the same V to the 
correlated site. 

For a cluster approximation with L c > 1, baths should 
be different for inequivalcnt correlated sites. It has to 
be expected, for example, that bath sites at the cluster 
boundary are more efficient to compensate for the dis- 
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regarded inter-cluster hopping processes than bath sites 
coupled to the cluster center. This can be studied using 
the reference system G in Fig. [1] which includes a cou- 
pling to a bath site at the edges (Vi) and at the central 
sites (V2) of a cluster with L c = 4 correlated sites. Fig. 
[5] shows the results of the according VCA calculation. 
Both the inner and the outer bath sites couple to the 
system with a finite hybridization and thereby lead to a 
decrease of the optimal SFT grand potential as compared 
to vanishing hybridization. For the outer ones, however, 
this binding-energy gain is about one order of magnitude 
higher. Also the optimal hybridization is much larger for 
the outer bath sites. 

This suggests to neglect the coupling of bath sites to 
the correlated sites at the center of the chain completely, 
i.e. to switch off V2 and to consider reference system H 
in Fig. [T] As can be seen in Fig. this represents an 
excellent approximation. Comparing the results for the 
reference systems G and H with each other by looking 
at the trend of fi(Vi) = f2[X(Vi)] for optimal V% and 
for V2 — 0, respectively, we find the optimal SFT grand 
potential to be only slightly higher and the optimal hy- 
bridization V± almost unchanged. 

The idea of attaching bath sites at the cluster edges 
only is pursued with the calculations shown in Fig. [5] 
We employ reference system H. For any cluster size from 
L c = 2 to L c = 8, it is found that edge bath sites couple 
to the system and decrease the SFT grand potential. For 
L c = 4 and [7 = 4 (see figure) this decrease amounts to 
Afl/L w 0.002. For stronger interactions the cluster ap- 
proximation generally tends to improve. Consequently, 
the optimal hybridization becomes smaller. The opti- 
mal grand potential at U = 8, for example, decreases 
by Afl/L w 0.0001 due to the bath sites which is one 
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ping. Arrows indicate minima. Dashed black line: BA result 
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order of magnitude less than for U = 4. It must be em- 
phasized that thanks to the Q-matrix technique^ which 
completely avoids frequency summations or integrations, 
there are no numerical problems to accurately compute 
energy differences of this or even lower order of magni- 
tude. From the numerical point of view, this is an im- 
portant step forward as compared to earlier evaluations 
of the SFT grand potential using integrations over fre- 
quencies to + irj with real to and subsequent extrapolation 
77 — >• (cf . Refs. I2DI34I for a detailed comparison) . 

Attaching bath sites and thereby allowing the electrons 
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to hop into the bath and back to the original site and 
thereby to gain kinetic energy, turns out to be more ef- 
fective than the gain in kinetic energy that is obtained by 
optimizing (and increasing) the intra-cluster hopping t' . 
This is demonstrated in the figure by calculations using 
the reference system I where the hybridization to edge 
bath sites and the intra-cluster hopping arc optimized 
simultaneously Comparing the results for reference sys- 
tems H and I shows that for any L c the binding energy 
gain due to the coupling of bath sites is considerably 
larger. 

Fig. [9] also includes the result of a calculation using 
reference system J (with L c = 4) where there is one bath 
site for any correlated site (not only at the edges) but still 
only a single hybridization strength that is optimized by 
assuming this to be the same for all bath sites. It is 
interesting to note that this reference system turns out 
to be inferior as compared to I and also to H (the SFT 
grand potential at the minimum is higher) although there 
are two more bath sites. This once more demonstrates 
the ineffectiveness of bath sites coupled to the center of 
the Hubbard chain. 

With increasing L c the optimal SFT grand potential 
(using reference system H, for example) nicely converges 
to the exact value which is shown in Fig. [5] for compar- 
ison. It is important to note that the inclusion of bath 
sites hardly speeds up this convergence. For any given 
cluster size L c , the additional inclusion of two more bath 
sites gives a binding-energy gain considerably smaller 
than the gain obtained by a cluster with two more cor- 
related sites. This also holds true if more bath sites are 
taken into account. The decisive lowering of the SFT 
grand potential is always due to a larger cluster size. 
Concluding, bath sites are quite ineffective as far as the 
grand potential or the ground-state energy is concerned. 

It is an interesting question whether or not bath sites 
at the edges of the cluster finally decouple from the corre- 
lated sites for L c — > oo, i.e. whether or not the respective 
optimal V vanishes in this limit. For bath sites coupled 
to the center of the cluster, a decoupling V — > is quite 
plausible physically and is actually foreshadowed by the 
results shown in Fig. [5] for reference system G. For edge 
bath sites, the optimal V is shown in Fig.llOlas a function 
of the inverse cluster size. As a simple cubic spline ex- 
trapolation shows, the results are compatible with a finite 
V for L c — > co (red dashed line) but also V — cannot 
be excluded. From the results shown in Fig. [5] it appears 
obvious that a minimum of f2[S(i')]/L is found for any 
finite L c . A finite position of the minimum for L c — > oo 
would imply an interesting behavior of the SFT grand po- 
tential as function of V as this must become completely 
flat (at least in a finite V range around V = 0). 

While bath sites are of minor importance as concerns 
static quantities such as the ground-state energy, they 
are decisive for dynamical quantities and for the single- 
particle excitation gap A in particular. This shall be 
demonstrated in the following. As argued above, a cluster 
method is likely to fail close to the critical point U = 0. 
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TABLE I: Single-particle excitation gap at U = 4 as obtained 
within the VCA using different reference systems as indicated 
(see Fig. [TJ. We set t = 1. C-DMFT result for L c = 2 and 
n s = 3 (Ref.HJ). 



Predicting the gap away from the critical point, however, 
can serve as any strong test for a cluster technique. 

Tab. H] shows results for the A at intermediate cou- 
pling U = 4 as obtained within different approaches. By 
evaluation of Eq. ([3]) we find the exact value A = 1.287. 
As it is well known, a metallic state with A = is pre- 
dicted by static mean-field theory and, for U = 4, by 
dynamical single-site approximations (dynamical impu- 
rity approximation^ and DMFT). On the other hand, 
even the most simple cluster approach, i.e. the VCA for 
L c = 2 and n s = 1 (without bath sites), predicts a finite 
gap but strongly overestimates its size. Improvement is 
possible using larger clusters, but even for L c = 10 the 
gap is overestimated by about 18% (see also Fig.[|5]). This 
remains essentially unchanged even with a full optimiza- 
tion of 5 different hopping parameters (Tab. HI VCA for 
reference system C). 

Introducing bath sites completely changes the situa- 
tion. Using reference system H, i.e. L c = 2 correlated 
sites with one bath site attached to each (n s = 2), yields 
a gap which is drastically too small. The result becomes 
even worse, namely the gap is almost closed, when in- 
creasing the cluster size to L c = 4 but still keeping one 
bath site attached to each of the two edge correlated sites. 
We infer that while the n s = 2 cluster nicely improves 
the ground-state energy it apparently fails to describe 
the excitation gap. 

Using one more bath site (reference system K, n s = 3) 
and L c = 2 yields a further but negligibly small decrease 
of the ground-state energy but a gap that comes very 
close to the exact one (which is underestimated by about 
8%). Now, an increase of the cluster size to L c = 4 yields 
further improvement, namely a gap that underestimates 
the exact one by 6% only. Adding another bath site at 
each of the cluster edges, i.e. n s = 4 (reference system 
L), yields an almost vanishing gap again. We conclude 
that there is a sizable odd/even effect with respect to the 
number of bath sites, and that a reliable prediction of the 
gap requires an even number (i.e. n s odd). Given this, 
the inclusion of bath sites is of crucial importance for an 
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accurate determination of single-particle excitations and 
the insulating gap. 

We would like to stress that due to the SFT variational 
principle and due the Q-matrix technique there is no ad- 
justable parameter in the calculation of the gap, once the 
reference system is specified. It has been verified that the 
results are converged with respect to the size of the orig- 
inal D = 1 Hubbard model (L ~ 10 3 ) and with respect 
to the number of Lanczos steps (Sl ~ 100). Within the 
cellular DMFT (using Lanczos as a cluster solver) on the 
other hand, the gap value somewhat depends on the pro- 
jection criterion employed to fix the bath parameters. In 
Ref. [3(3 the definition of the C-DMFT gap is furthermore 
adjusted to recover the exact value in the strong-coupling 
limit. Thereby, a gap of A = 1.14 is obtained for an 
L c = 2, n s = 3 cluster which is close to our result (see 
Tab.rj|. 

Using the reference system K with L c = 4 and n s = 3, 
we have computed the gap also for different U. The de- 
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viation to the exact result is found to decrease further 
for stronger interactions and amounts to 1.7% for U = 6 
and 0.3% for U = 8. The gap size is underestimated. For 
smaller U, deviations are larger (underestimation of 10% 
and 7% for U = 3 and (7 = 2, respectively). Compared 
with the C-DMFT results^ however, this is still a con- 
siderable improvement. As a cluster approach the VCA 
cannot reproduce the exponentially small gap in the crit- 
ical regime and finally overestimates the gap size with a 
relative error of about 340% at U = 1 being a typical 
example (here the exact gap is A = 0.005). 

To analyze the origin of the apparent odd/even ef- 
fect, we discuss the interacting local density of states 
(LDOS) which is shown in Fig.[TT]for (7 = 4. The LDOS 
is calculated from the imaginary part of the local VCA 
Green's function which, opposed to DMFT or C-DMFT, 
differs from the cluster Green's function. Since a real- 
space cluster approach necessarily breaks translational 
symmetry, the LDOS is computed for a central cluster 
site where convergence for L c — > oo is expected to be the 
fastest. The LDOS turns out to be non- vanishing on a 
large number of small but finite frequency intervals. This 
structure is an artifact that is due to the finite Hilbcrt 
space corresponding to the reference system and due to 
the corresponding discrete pole structure of the VCA self- 
energy. A smooth LDOS can therefore only be obtained 
with some additional broadening. For the figure we have 
used a Lorentzian broadening with a comparatively large 
broadening parameter rj = 0.2 (blue lines) while the gap 
size (which is calculated for 77 = 0) can be read off quite 
accurately from the spectra with small 77 = 0.005 (red 
lines) . 

Comparing the results for the different reference sys- 
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terns with each other, we find the same overall structure 
of the LDOS in all cases but also significant finite-size 
effects. The latter prevent a reliable prediction of the de- 
tailed shape of the LDOS. While we expect the L c = 10 
cluster (without bath sites, upper panel) to give the best 
estimate for the LDOS in general, the discussion above 
has shown that the L c — 4, n s — 3 cluster yields a much 
more reliable prediction of the insulating gap (middle). 
Looking at the result for -q = 0.005, a deviation from the 
exact result is in fact hardly visible on the scale of the 
figure for this cluster. Contrary, the LDOS computed for 
the L c = 4, n s = 2 cluster (lower panel) clearly shows 
finite spectral weight for frequencies much smaller than 
the exact gap (see also the inset). This weight is rather 
small but significant as can be seen by varying the chem- 
ical potential: A slight decrease from its particle-hole 
symmetric value [i = J7/2 of about Afi = 0.04 is suffi- 
cient to produce a metallic state. 

To understand the failure of the L c = 4, n s = 2 ap- 
proximation, one can (ad hoc) decrease the intra-cluster 
nearest-neighbor hopping from t' = t to t' = in the 
reference system. This results in an LDOS with a three- 
peak structure consisting of the two Hubbard bands and 
a quasi-particle peak at to = as it is typical for the 
metallic solution of the half-filled Hubbard model within 
a dynamical impurity (mean-field) approximation given 
by L c = 1. The spectral weight in Fig. [11] (lower panel) 
that is responsible for the too small gap in the case of 
n s = 2 (odd number of bath sites) , thereby continuously 
evolves into the quasi-particle peak of the metallic mean- 
field LDOS. Now, within the dynamical impurity approx- 
imation, such a metallic LDOS at half-filling is always 
produced by an odd number of bath sites where one bath 
site has a one-particle energy exactly at the Fermi edge /i 
while the energies of the remaining are located symmetri- 
cally to fi. Due to a decreasing and eventually vanishing 
hybridization strength, this Fermi-edge bath site decou- 
ples from the rest of the reference system upon approach- 
ing the Mott insulator with increasing U. The Mott insu- 
lator itself is therefore characterized by an even number 
of bath sites. We therefore interpret the occurrence (ab- 
sence) of low-frequency spectral weight within the exact 
gap as being a reminiscence of the low-frequency charac- 
ter of the corresponding metallic (insulating) mean-field 
solution. Hence, a reasonable description of the Mott 
insulator and a reliable prediction of the insulating gap 
requires an even number of bath sites per correlated site 
also in the case of cluster (L c > 1) approximations. 

Finally, we compare our results for the local Green's 
function on the imaginary frequency axis with calcu- 
lations by Hafermann et al. (see Ref. l45h for U = 6. 
Fig. [T2] shows our VCA results for the L c = 2 cluster 
with n s = 3 local degrees of freedom per correlated site 
in comparison with essentially exact dynamical DMRG 
data. As compared with L c = 2 cellular DMFT using 
weak-coupling continuous-time QMC at finite tempera- 
ture (1/T 1 = P = 20), there is a marginal improvement 
only. Wc would like to stress, however, that the VCA cal- 



culations are computationally much cheaper. The cluster 
dual-fcrmion approach for L c = 2 (and j3 = 20) gives a 
considerably better result and is already quite close to 
the DMRG data. Further improvement is possible for 
larger clusters. VCA calculations for the L c = 4, n s = 3 
reference system, i.e. with two more correlated sites, are 
hardly distinguishable from the DMRG. 



VI. FILLING-DEPENDENT MOTT 
TRANSITION 

To complete the benchmarking of different cluster ap- 
proximations, we study the metallic phase off half-filling 
and the filling-dependent Mott transition. The metallic 
phase is characterized by a finite compressibility k (see 
Eq. (H))). Changing the electron density (filling) n by 
changing the chemical potential fi at fixed U, the Mott 
insulator at n — 1 is approached with a diverging n — > oo 
for n -> 1 (see Eqs. © and ©). 

An approximation that is based on a cluster with a 
finite (and small) number of degrees of freedom, neces- 
sarily implies a strongly limited frequency resolution and 
thus a rough description of low-frequency physical prop- 
erties. We therefore expect that it is more difficult within 
the VCA to describe a metallic state with gapless single- 
particle excitations as compared to the description of a 
gapped Mott-insulating state. On the other hand, this 
argument neglects the fact that the limitation of the fre- 
quency resolution can at least partly be compensated for 
by an adaption of the cluster one-particle parameters. 
This is well known from the DMFT-ED approach^ as 
well as from dynamical impurity approximations^ which 
show that the low-frequency quasi-particle (Kondo) res- 
onance can be accessed with a few sites only. At least 
technically, however, VCA calculations arc more difficult 
for the metallic state. This is simply due to the fact 
that, for a given cluster size, the absence of particle-hole 
symmetry implies an increased number of independent 
variational parameters to be optimized. 



A. Variation of one-particle energies 

It suggests itself, for example, to include an overall 
shift of the one-particle energies of the cluster sites into 
the set of variational parameters. Note that for the 
particle-hole symmetric case this shift is irrelevant (the 
VCA grand potential is stationary at a vanishing shift). 
At a finite hole concentration, however, the grand po- 
tential is at a maximum for a finite shift (which is dif- 
ferent from the shift of the chemical potential) . As has 
been shown in Ref. l34l . this ensures the thermodynamical 
consistency of the approach with respect to the particle 
number. Especially for the filling-dependent Mott tran- 
sition, it is important that the filling calculated from the 
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FIG. 13: (Color online) Electron density n c = {N)/L c = 
N c /L c as a function of the chemical potential for an L c — 10 
Hubbard chain (open boundary conditions) compared to the 
exact result for L c — > oo from Ref. [23. 



(approximate) interacting density of states, 
1 /"° 

« = j ^2 / doj p la (uj) , (12) 

where /Oi CT (w) = (— l/7r)ImG iicr (tJ + i0 + ), gives the same 
result as the filling calculated from the (approximate) 
grand potential, 



1 dfl 
L dp 



(13) 



i.e. that calculations on the level of one-particle exci- 
tations are consistent with those on the (zero-particle) 
static thermodynamical level. 

The one-particle excitation spectrum is most accu- 
rately determined by using the dynamical density-matrix 
renormalization-group (DMRG) technique^ or the quan- 
tum Monte-Carlo approach. 47 i 48 The essential features 
of the spectrum, including the more intense spinon and 
holon bands, however, are already accessible using fi- 
nite Hubbard chains of rather small size (e.g. L c = 10). 
This has already been demonstrated by calculations us- 
ing cluster-perturbation theory (CPT)^ i.e. VCA with- 
out any parameter optimization at all. Even for L c = 2 
and a number of n s — 1 = 2 additional bath sites per cor- 
related site, the overall spectrum is in very good agree- 
ment with the more accurate DMRG results as has been 
demonstrated by Capone et ali^I 

Here we would like to focus on a different point which 
is relevant for any cluster approach. This is illustrated in 
Fig. [T3] for an isolated cluster with L c = 10 sites. Due to 
the U(l) symmetry of the grand-canonical cluster Hamil- 
tonian H', the cluster ground state has a definite total 
particle number N c which, depending on the chemical po- 
tential /i, can vary from iV c = to N c = 2L C . Hence, the 
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FIG. 14: (Color online) Electron density n within cluster- 
perturbation theory (CPT), calculated by integrating the 
CPT density of states up to zero excitation energy (red line, 
Eq. (|12p ) and calculated by differentiation of the CPT grand 
potential with respect to fi (blue line, Eq. (|13(l ). as function of 
(j, compared to the exact result (BA). Cluster size: L c — 10. 
Horizontal dashed lines: cluster fillings. Vertical dashed lines: 
critical chemical potentials at which the cluster ground state 
changes. 



density n c = N c /L c acquires discrete values only when 
varying fj. and discontinuously jumps at certain critical 
chemical potentials fi c j. As is demonstrated by the fig- 
ure, this is a strong finite-size effect that cannot be tol- 
erated if one is interested in the filling-dependent Mott 
transition of the system in the thermodynamical limit. 

Therefore, it is an obvious question whether it is possi- 
ble to predict a continuous and reliable trend of the filling 
as a function of /z for L c — * oo but using an approxima- 
tion based on finite (small) clusters only. This is also 
related to the question whether one can access systems 
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FIG. 15: (Color online) The original system and different 
reference systems with corresponding variational parameters. 
See text for discussion. 
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with an arbitrarily low hole concentration S = 1 — n as it 
is necessary, for example, to recover the compressibility 
divergence for 5 — > 0. 

By considering an infinite system of disconnected clus- 
ters and by re-introducing the inter-cluster hopping in 
the lowest non-trivial order, cluster-perturbation theory 
directly works in the thermodynamical limit. Therefore, 
CPT could by able to predict a continuous \i dependence, 
in principle. As can be seen from Fig. [T4] CPT in fact 
gives a metallic state with a compressibility k that is fi- 
nite everywhere except for the critical chemical potentials 
of the isolated cluster at (i = fj, Ct i. This is a substantial 
improvement as compared to the direct diagonalization 
where K = (except for [i = (J, c ,i)- However, the obvi- 
ous disadvantage is that the CPT filling exhibits finite 
jumps at [i = fj, Ci i. This is easily understood by looking 
at Eq. (|12p since for each /i = fi c ^ there is a discontinuous 
change of the ground state of the (grand-canonical) clus- 
ter Hamiltonian which implies a discontinuous change of 
the cluster Green's function G'(ui) and thus of the lattice 
(CPT) Green's function G(w) = {G'^y 1 - (t - t'))' 1 . 

Essentially the same applies to the case where n is 
calculated as the \i derivative of the grand potential, see 
Eqs. (fT0|) and (fT3|) . Furthermore, the comparison of the 
results obtained from Eqs. (fT2|) and (fl"3]l illustrates the 
thermodynamical inconsistency of the CPT (see Fig. [T4"]) . 
Compared with the exact Bethe ansatz result, it turns 
out to be much better to calculate the filling from the 
CPT density of states. 

We conclude that plain CPT cannot describe the 
filling-dependent Mott transition and exhibits severe 
problems in describing the trend of n as a function of 
the chemical potential. Restricting the approach to the 
discrete set of cluster densities, however, gives rather sat- 
isfactory results. Fig. PHI shows that the chemical poten- 
tial, where the CPT filling equals one of the accessible 
cluster fillings, is close to the exact /i corresponding to 
that filling and that both ways to compute n (Eqs. fT2]) 
and (fT3]) ) almost yield the same result. This is nicely 
demonstrated in the inset for n = 0.4 and n = 0.6. How- 
ever, there are still problems, even for an accessible clus- 
ter filling as, for example, in the case of n = 0.5 where 
no unique chemical potential can be read off and n = 
in a finite fi range (see inset again). 

The situation changes completely when using the vari- 
ational cluster approximation. The most simple reference 
system consists of isolated clusters where only an overall 
shift of the on-site energies is taken as a variational pa- 
rameter to guarantee thermodynamical consistency, see 
Fig. fTS] A. The corresponding results for clusters of dif- 
ferent size L c are shown in Fig. [16] 

Apart from thermodynamical consistency, the most 
important difference as compared to the CPT, consists 
in the fact that a continuous /i dependence is found in 
the entire range from fx = — 2 (empty band) to fi = fj, c 
(half-filling). This has explicitly been checked for clus- 
ters with L c = 2, 4, 6, 8, 10. It turns out that the optimal 
value for the shift of the on-site energies (partly) com- 
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FIG. 16: (Color online) Electron density (filling n) as a func- 
tion of the chemical potential close to half-filling for U = 4 as 
obtained from the VCA with reference system A of Fig. ll5l for 
different cluster sizes L c . In all cases the cluster ground state 
is in the N c = L c subspace. Inset: n{fj) for L c = 4 displayed 
for the entire /j, range and with cluster ground state in the 
N c — L c and in the N c — L c /2 subspace. 



pensates for the deviation of the chemical potential from 
its particle-hole symmetric value /j, = U/2 such that (due 
to the large Mott-Hubbard gap) the cluster ground state 
is always found in the subspace where the cluster itself 
is half-filled (N c = L c ). The cluster ground state, the 
optimal self-energy and eventually the filling n thereby 
continuously depend on (i. This mechanism also works 
down to fi = —2. 

Fig. HU displays the critical regime close to half-filling 
only. We find that the critical value of the chemical po- 
tential fj, c where the transition from the metal to the 
Mott-insulating state takes place, increases with increas- 
ing cluster size L c and converges to the exact result. The 
figure also shows, however, that this convergence is rather 
slow. Furthermore, it remains unclear whether or not the 
compressibility divergence can be recovered in the limit 
L c — * oo: Note that the slope of n(/n) (i.e. k) for n — ► 1 
appears to decrease with increasing L c . This might be 
explained by the argument that, even for L c — > oo and 
even with optimized cluster one-particle parameters, a 
trial self-energy taken from the N c = L c subspace cannot 
describe the physics of a metallic state entirely correct. 

We conclude that the continuous dependence on /i is 
achieved at the cost of fixing the cluster ground state 
to half-filling. This problem becomes more and more 
severe with decreasing filling. The inset of Fig. [16] shows 
a calculation for L c — 4 illustrating this issue. For the 
calculation with 7V C = L c = 4, the filling is reasonably 
close to the exact filling in the vicinity of half-filling only. 
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In the vicinity of quarter filling, however, a much better 
result is obtained with a VCA calculation starting from a 
cluster ground state with N c = 2 = L c /2 (quarter- filled 
cluster). This is physically plausible. 



B. Bath degrees of freedom 

As has been seen in the discussion of the results for 
half-filling, bath sites can considerably help to improve 
a cluster approximation. This is all the more important 
in the case of a metallic system off half-filling since bath 
sites also serve as charge reservoirs. Varying the chemical 
potential or another physical model parameter, the elec- 
tron density on the correlated sites can vary smoothly by 
a charge flow from and to the uncorrelated bath sites in 
the reference system. 

Opposed to (cellular) DMFT, the filling in the origi- 
nal model n, as calculated from the approximate VCA 
lattice Green's function, is generally different from the 
density at the correlated sites in the (cluster) reference 
system n' . By rule of thumb, however, the deviations 
are small, i.e. n«n'. This implies that all fillings from 
n = up to n = 1 (half-filling) can be realized by using 
a strictly half-filled reference system provided that the 
cluster includes (at least) one bath site per correlated site 
(n s = 2). Consider, for example, a cluster with L c = 2 
correlated and n s — 1 = 1 bath site per correlated site. 
In total the cluster then consists of L' = 4 sites. For 
particle-hole symmetric parameters, the cluster ground 
state lies in the subspace with N(. ot = 4 electrons, and 
symmetry arguments imply an electron density n' = 1 
on the correlated and n' hath = 1 on the bath sites. This 
corresponds to half- filling, n = 1, for the original model. 
For n < 1 wc will find n' ~ n < 1 and n' heith > 1 such that 
n '+ n bath = 2, i.e. a half-filled cluster ground state. In the 
limit n — > 0, the iV t ' ot = 4 electrons will mostly be located 
on the bath sites, i.e. n — > and n' hllth — > 2. Analogous 
arguments hold for fillings above half-filling. This mech- 
anism promises continuous dependencies on the chemi- 
cal potential with a cluster ground state staying in the 
iV t ' ot = V subspace, while the physical properties are 
governed by the density on the correlated sites n' ~ n, 
which varies smoothly with /i. 

We have not been able to find a stationary point by us- 
ing two bath sites per correlated site (n s = 3). Remem- 
bering the discussion of the half-filled case, this does not 
appear to be uncommon, since a reference system with an 
even number of bath sites per correlated site is expected 
to give a good description of the Mott insulator but not 
of the metal. For the calculations we therefore concen- 
trate on reference systems with n s = 2. For simplicity, 
we attach one bath site to any of the correlated sites and 
assume the hybridization V to be the same for all sites. 
Additional variational parameters are e and e' , the one- 
particle energies of the correlated and of the bath sites 
which arc assumed to be constant again. The reference 
system is displayed in Fig. [T5l B. 
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FIG. 17: (Color online) Filling n as a function of [i close to 
half-filling for U = 4 as obtained from the VCA with reference 
system B of Fig. [15] for L c = 2 and L c = 4. Exact and C- 
DMFT result from Ref. [3l] are shown for comparison. 



Numerical results for clusters with L c = 2 and L c = 
4 correlated sites are shown in Fig. [T7l Irrespective of 
the cluster size, there is an excellent agreement with the 
exact result for fillings lower than n ~ 0.75. But also 
for higher fillings the VCA results with bath sites are 
convincing and represent a considerable improvement as 
compared to the CPT results (Fig. IT¥|) but also compared 
to the VCA results without bath sites (Fig. [TBJ>. The 
critical chemical potential [i c for the transition to the 
Mott insulator is somewhat overestimated but the error is 
much smaller than the underestimation of fj, c within VCA 
without bath sites. More important, however, it appears 
that the approach correctly predicts the divergence of the 
compressibility. Unfortunately, it has turned out to be 
extremely difficult numerically to follow up the stationary 
point as a function of fi in the region very close to half- 
filling. The C-DMFT results of Capone et ali^i which arc 
also shown in the figure, arc slightly closer to the exact 
n(/i). Note, however, that this has been achieved and 
crucially depends on a special (but physically motivated) 
choice for the distance measure which emphasizes the low 
Matsubara frequencies. In contrast, our approximation 
is free from any adjustable parameter. 

Note that L c and the number of bath sites per corre- 
lated site (i.e. n s — 1) do not fully specify the reference 
system. Different systems can be generated by the differ- 
ent ways in which bath sites are coupled to the correlated 
ones. Reference system C in Fig. 1151 for example, is char- 
acterized by L c = 2 and n s = 2 but spans (when inde- 
pendently varying all one-particle parameters) a space of 
trial self-energies which differs from the one spanned by 
reference system B (with L c = 2). In the limit n s — ► oo, 
i.e. for the case of cellular DMFT, the different ways of 
coupling the baths to the correlated sites do not matter 
as they can be mapped onto each other by unitary trans- 
formations and therefore span the same space of trial 
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FIG. 18: (Color online) Chemical-potential dependence of the 
optimal one-particle parameters for reference system C of Fig. 
1 1 5 1 at U — 4. Ae is an overall shift of all one-particle energies 
(e, £ a and e b ). 
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FIG. 19: (Color online) Filling versus chemical potential for 
U — 4. VCA result using reference system C of Fig. [15] in 
comparison with the exact result (Bethe ansatz), with DMFT 
(from Ref. W0 ) and with the two-site dynamical impurity ap- 
proximation (DIA). 



sclf-encrgics. This is different, however, for small n s . As 
is demonstrated in the following, reference systems B and 
C yield very similar results for the metallic phase while 
C gives a much better description of the Mott insulator. 

Any bath site in reference system C is connected via 
hybridizations V\ and V2 to both correlated sites. Re- 
quiring the self-energy to be symmetric with respect to 
an interchange of the two correlated sites, implies that 
the modulus of the two hybridization parameters must 
be the same, i.e. V\ = 1LV2 (if the hybridization is as- 
sumed to be real). We consider two bath sites, one with 
V\ = — V2 = V a and another one with V\ = V2 = V&. 
This is the only choice left if the reference system is re- 
quired to respect particle-hole symmetry for /i = U/2. 
Consequently, there are five independent variational pa- 
rameters in total, V a and Vb, the bath on-site energies e a 
and Eb and a general shift of all on-site energies (including 
the correlated sites) Ae. 

Fig. [18] shows the optimal values of these parameters 
as functions of the chemical potential. At half-filling in 
the Mott insulator, i.e. for \x > \i c , the parameters are 
fi independent. The overall shift of the on-site energies 
Ae vanishes, V a = Vb, and e Q + £b ~ U. For the particle- 
hole symmetric point at fi = U/2, these restrictions are 
enforced by symmetry and respected by the approxima- 
tion. The parameters evolve continuously with \x except 



for /i = /i c « 1.42, i.e. at the metal-insulator transition. 
With decreasing [i, the bath on-site energies decrease. 
This results in an increasing bath-site density n' and thus 
in a decreasing filling n as described above. In the en- 
tire \x range, the reference system remains strictly half- 
filled. At 11 ?s —0.76 the hybridization of bath site "a" 
vanishes, i.e. this site decouples from the system. Since 
E a + Ae < fj, the site is completely occupied. For lower /j, 
the remaining reference system therefore consists of the 
two correlated and merely a single uncorrelated bath site 
to be filled with two electrons in total. As n — * —2 also 
the second bath site "b" decouples {Vb — ► 0) taking both 
electrons left such that the filling n — > 0. The bath sites 
perfectly do their job as charge reservoirs. 

The resulting filling as a function of /i is displayed in 
Fig. [TO] The agreement with the exact Bethe ansatz re- 
sult is excellent for low and intermediate fillings. A slight 
deviation is found for fillings n > 0.8. Nevertheless, the 
qualitative trend is predicted correctly including the di- 
vergence of the compressibility. As compared to the re- 
sults of reference system B, there is hardly any difference 
on the scale of the figure. 

The critical chemical potential turns out to be fi c ~ 
1.42 within the VCA. This should be compared with the 
VCA result (using the same reference system) for the 
Mott-insulating gap A = 1.128 which is read off from 
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the LDOS at half-filling. Via 

A = 2(t//2 - Mc ) (14) 

this implies a critical chemical potential of fi c = 1.436. 
The small difference in the results for \i c is most probably 
caused by numerical problems to locate the stationary 
point for the critical regime in the metallic phase near 
the transition. The predicted gap is close to the gap 
found with reference system K in Fig. [T] for L c = 2 which 
implies that, opposed to reference system B (Fig. [15]), 
not only the metallic phase in the entire filling range but 
also the Mott insulator is well described. Considering 
the small size of the cluster this is actually a surprisingly 
good and very satisfying result. 

We also note that the VCA calculation fully respects 
Luttinger's sum rule*^ in the entire filling range. While 
for a translationally invariant system the sum rule states 
that the filling equals the range in reciprocal space en- 
closed by the Fermi points, a generalized form of the 
sum rule (given in Ref. 1351 ) must be considered here 
to account for the artificially reduced symmetry result- 
ing from the cluster approximation. The main idea is 
that the sum rule can be derived from the equation 
liniT^o Tr(S<9G / d(iu) n )) = 0. This can be tested even if 
the approximate self-energy and the (via Dyson's equa- 
tion) related Green's function do not respect the trans- 
lational symmetries of the underlying lattice and even if 
the system under consideration is finite. We find that the 
sum rule already holds for the reference system C itself. 
As has been discussed recently^ due to the conserving 
nature of the VCA, the validity of the sum rule for the 
reference system is then transferred to the approximate 
quantities for the lattice problem. 

Fig. [TO] also includes the prediction of single-site 
DMFT. While close to half-filling the mean-field theory 
fails completely, hardly any difference to the exact re- 
sult can be detected for fillings lower than n « 0.6. It 
is worth mentioning that the results of full DMFT arc 
quantitatively recovered by the most simple dynamical 
impurity approximation (DIA) within the framework of 
SFT, namely by using a reference system consisting of 
the single correlated and a single bath site only (L c = 1 , 
n s = 2). Both n(fi) curves, taken from DMFT and from 
the DIA, arc identical on the scale of the figure. 

VII. CONCLUSIONS 

Cluster mean-field theories that arc based on full diag- 
onalization or on the Lanczos technique to treat the ef- 
fective cluster problem, exhibit a number of advantages: 
They directly work at zero temperature, they are flexi- 
ble and can equally well treat arbitrary geometries, and 
they provide the numerically exact solution of the cluster 
within a comparatively short CPU time. These advan- 
tages are achieved at the cost of a strongly limited cluster 
size (of the order of 10 sites) dictated by the exponen- 
tial dependence of the cluster Hilbcrt-spacc dimension on 



system size. It is therefore of highest importance to make 
use of the cluster degrees of freedom in the best possi- 
ble way. This can be accomplished with the self-energy- 
functional theory. The SFT allows to fix the cluster one- 
particle parameters with the help of a physical variational 
principle that is constructed for the optimization of the 
one-particle excitation properties. 

The actual choice of the reference system, however, is 
not prescribed by the approach itself, i.e. different clus- 
ter topologies and thus different cluster approximations 
are conceivable. This is the motivation for the present 
study. With the focus on the interaction- and filling- 
dependent Mott metal-insulator transition, as a prime 
example of a correlation effect, different reference sys- 
tems have been tested against each other and against ex- 
act results available for the one-dimensional case. Note 
that a one-dimensional lattice model actually represents 
the most difficult test case for a cluster approximation. 

In the following we recapitulate the main results of our 
study. First of all, the Mott insulating state of the model 
at the particle-hole symmetric point is well described by 
a cluster approximation using a few correlated sites only. 
In particular, the ground-state energy can be determined 
precisely using finite-size scaling. Already the most sim- 
ple implementation of the variational-cluster approach 
(VCA) including merely an overall optimization of the 
intra-cluster hopping, clearly improves on the "direct" 
cluster approximation. However, as compared to plain 
cluster-perturbation theory without any parameter op- 
timization at all (but combined with the expression for 
the SFT grand potential), the VCA yields a marginally 
improved ground-state energy only. 

Similarly, the independent optimization of several or 
even of all intra-cluster nearest-neighbor hopping param- 
eters leads to a gain in binding energy but this is small 
compared to the gain obtained by increasing the clus- 
ter size by two more sites. With increasing distance 
to the chain edges, the optimized hopping quickly con- 
verges to the "physical" value, and bulk properties are 
already found for sites at a distance of more than two 
nearest-neighbor units from the cluster edge. Hopping 
parameters vanishing in the original model can acquire a 
non-zero but small value in the optimized reference sys- 
tem. Hopping parameters linking the edges of a cluster 
arc found to vanish, i.e. the VCA prefers open bound- 
ary conditions. The same holds for hopping parameters 
which would imply a breaking of particle-hole symmetry: 
A next-nearcst-ncighbor hopping, for example, turns out 
to be zero at the stationary point. 

The VCA correctly predicts a Mott insulating state for 
any U > 0. It also gives a reasonable estimate for the 
size of the single-particle insulating gap. This estimate 
improves with increasing interaction strength. While for 
strong and intermediate coupling the optimal results (for 
the ground-state energy as well as for the gap) are ob- 
tained for reference systems with parameters close to the 
original model, strongly deviating parameters (more than 
100%) are favorable in the weak-coupling regime. Even- 
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tually, for U — > the VCA (and presumably any real- 
space cluster approach) fails to describe the low-energy 
physics of the Mott transition. Remembering the cluster 
mean-field nature of the approximation, this had to be 
expected. It is obvious that critical behavior cannot be 
accessed while, on the other hand, it is satisfying that 
physical properties on a higher energy scale are accu- 
rately captured with rather small clusters only. 

The additional consideration of bath sites in the refer- 
ence system always yields an improved ground-state en- 
ergy. It has turned out that bath sites tend to decouple 
from a cluster reference system at the correlated sites in 
the cluster center while they tightly couple to the system 
at the chain edges. As concerns the ground-state energy, 
however, this hardly speeds up the convergence to the 
exact result with increasing cluster size. Going to the 
next larger cluster is always found to be more effective. 

On the contrary, bath sites decisively influence the de- 
scription of the one-particle excitation spectrum. Cluster 
reference systems with an even number of additional bath 
sites coupled to a correlated site (at the cluster edge) give 
a much better result for the gap than clusters without 
bath sites but more correlated sites. This demonstrates 
the importance of local correlations for the one-particle 
spectrum. An odd number of bath sites per correlated 
site overemphasizes the mean-field character and thereby 
fails to predict the gap accurately. 

The overall features of the local density of states (in- 
sulating gap, moments, etc.) are addressed when looking 
at the local one-particle Green's function on the imag- 
inary frequency axis. Using essentially exact dynami- 
cal DMRG results as a benchmark, we find our VCA 
results to be fully competitive with cellular DMFT cal- 
culations for the same L c . While C-DMFT combined 
with a stochastic technique as a cluster solver is much 
more time-consuming, C-DMFT (L c = 2, n s = 3) com- 
bined with Lanczos is less reliable as the VCA (L c = 2, 
n s = 3). The latter gives a significantly better result for 
the insulating gap which, apart from the choice of the 
reference system, is unbiased and free of parameter fit- 
ting. This shows that the thcrmodynamically consistent 
determination of the cluster parameters is worth the ef- 
fort. Predicting the detailed shape of the local density of 
states (for real frequencies) is beyond the capabilities of 
a cluster approach based on the Lanczos method: Finite- 
size effects are clearly present . 

As concerns characteristic quantities of the Mott insu- 
lator at half-filling, one can state that the convergence 
to the exact result is mainly determined by the number 
of correlated sites. The additional inclusion of a contin- 
uum of bath degrees of freedom or of a large number of 
bath sites appears to be unnecessary. A few bath sites, 
however, can strongly improve the results. For the de- 
scription of the metallic phase off half-filling and of the 
filling-dependent Mott transition this is even more cor- 
rect: 



Without the inclusion of bath degrees of freedom, the 
VCA with optimization of an overall on-site energy shift 
predicts a smooth dependence of the filling n on the 
chemical potential /j,. While this represents a clear ad- 
vantage as compared to the "direct" cluster approach and 
to cluster-perturbation theory, the main problem consists 
in the fact that the cluster particle number remains con- 
stant. Consequently, starting from a half-filled cluster, 
results are less and less reliable for increasing hole con- 
centration 1 — n. 

On the other hand, for reference systems with at least 
a single bath site per correlated site (n s > 2) there is 
a mechanism which solves the problem: With the total 
cluster particle number always being equal to the total 
number of cluster sites, iVtot = L', the particle density at 
the correlated sites and, related to that, also the filling 
of the original system may vary continuously in its entire 
range. 

For clusters with L c correlated sites and n s — 1 = 1 ad- 
ditional bath site attached to each, this mechanism has 
turned out to work over a wide range of fillings. Ex- 
cellent agreement with the exact n(/i) curve from the 
Bethe ansatz is found for fillings lower than n w 0.75. 
For fillings close to half-filling in the critical regime, the 
VCA still gives a qualitatively satisfying result, compara- 
ble to C-DMFT calculations, but slightly overestimates 
the critical chemical potential. While accessing the criti- 
cal regime for n — ► 1 poses difficulties, the compressibility 
divergence is clearly visible. 

A smooth \i dependence over the complete filling range 
as well as a good description of the Mott insulator can be 
achieved when using (for L c = 2 and n s = 2) a somewhat 
different cluster topology where a bath site couples to 
both correlated sites. We find that both the exact result 
as well as the best (DMFT) mean-field result are almost 
quantitatively recovered with small (L c = 2, L c = 1) 
reference systems including the minimum number (n s = 
2) of local degrees of freedom. 

The one-dimensional Hubbard model considered in 
this study has merely served for benchmarking pur- 
poses. Eventually, our main interest is focused on the 
physics of strongly correlated electrons on two- or higher- 
dimensional lattices. Doped Mott insulators in two di- 
mensions, however, constitute notoriously difficult many- 
body problems which require, as a prerequisite, thorough 
studies of systems which are more controlled in order to 
avoid artifacts or misinterpretations. 
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